A new representation of the bulk current in the quantum Hall effect regime 
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In preceding papers a Landauer-Biittiker type representation of bulk current transport has been 
successfully used for the numerical simulation of the magneto transport of 2-dimensional electron 
systems in the high magnetic field regime. In this paper it is demonstrated, that this representation 
is in full agreement with a treatment of the bulk current transport as a tunneling process between 
magnetic bound states. Additionally we find a correspondence between our network representation 
and the bulk current picture in terms of mixed phases mapped on a checkerboard: At half filled 
Landau level (LL) coupled droplets of a quantum Hall (QH) liquid phase and coupled droplets of 
an insulator phase phase exist at the same time, with each of them occupying half of the sample 
area. Removing a single electron from to such a QH liquid droplet at half filling completes the QH 
transition to the next higher QH plateau. Adding a single electron to such a droplet at half filling 
completes the QH transition to the previous lower QH plateau. As a consequence, the sharpness of 
the QH plateau transitions on the magnetic field axis depends on the typical size of the droplets, 
which can be understood as a measure of the disorder in the sample. 

PACS numbers: 73.43.-f, 73.40.Gk, 73.43. Nq, 73.43. Qt 
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I. INTRODUCTION 



Even more than 20 years after the discovery of the 
integer quantum Hall effect there are still controversial 
discussions about the origin of this phenomenon. In par- 
ticular the transition regime between plateaus gained in- 
creasing attention during the last two decades. While 
it is commonly accepted, that the plateau values of the 
integer quantum Hall effect (IQHE) can be explained by 
the edge channel (EC) picture 0, 0, 0, 0, the transition 
regime between plateaus is believed to be driven by a 
different mechanism in the bulk of the sample. A well 
defined scaling behavior of the temperature dependence 
is seen as an evidence for a quantum phase transition in 
the vicinity of the QHE plateau transitions 0, IE S H @ 
(for a review see also Refs. 0,^3). While the existence 
of a formal equivalence of the edge channel picture and 
the bulk current picture is proposed by Ruzin et al.[l^. 
the investigations of non-equilibrium situations like se- 
lective edge channel population, AC-characteristics and 
non-linear transport give evidence, that edge channels 
should be more than just a different representation of 
bulk currents. Since more than two decades the funda- 
mental question concerning edge vs. bulk in the QHE 
gets permanent attention even until now, as can be seen 
from a representative selection of references: Starting 
with the discovery of the IQHE ^| and first system- 
atic experiments investigating different contact configu- 
rations 0, 0] , continuing with the introduction of edge 
channels by Buttiker [l| and the enormous number of ex- 
amples for dealing with the edge vs. bulk problem and 
making use of the EC picture for mode ling and inter- 
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very fundamental aspect of edge channels is the drastic 
enhancement of the phase coherence length of electrons 
within edge states. Major work on this topic has been 
done by Komiyama et al [93| . 

Concerning modeling of quantum transport, one o f the 
very first attempts can be attributed to R.Landauer |ll4 
Ill2j |. A major step forward was achieved by M. Buttiker 
by introducing the so-called Landauer-Biittiker (LB) for- 
malism and the EC-picture of the QHE 1] . For modeling 
quantum transport and localzation in the bulk, mainly 
network model s on the basis of the Chalker Coddington 
(CC) network [ill are used. However, a model, which 
generates data in terms of voltages and resistances, which 
would allow a direct comparison which experimental data 
for realistically shaped samples had not been developed 
so far. In order to do so, an approach to combine EC 
transport and bulk transport has been made some time 
ago HE 0|- Subsequently that model has been expanded 
to a network [8^, which is also the subject of this paper. 
Although the basic idea of our network is common with 
the basic idea of the CC network, our handling of the 
nodes is substantially different: In contrast to a CC net- 
work our network does not use a transfer matrix for am- 
plitudes and phases. We use transmission by tunneling, 
but incorporating the effect of tunneling according to the 
LB formalism. The nodes are described by a back scat- 
tering function P, which is the ration of reflection and 
transmission coefficient R/T . 
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FIG. 1: Node of the network with two incoming and two out- 
going channels. The channels 1 — > 2 and 3 — > 4 are treated 
like edge channels with back scattering, where P = R/T ac- 
cording to the Landauer- Biittiker formalism. 




FIG. 2: Arrangement of the nodes for building the minimal 
physical element of a network, which is the closed loop of a 
so called magnetic bound state (see below) 



In FigQ] a single node of the network is shown, which 
transmits potentials from the incoming to the outgoing 
channels. The transmitted potentials of the out goin g 
channels as a function of P is calculated as follows 851: 



= (jM + P-pn)/(l + P) 



(la) 



(lb) 



In Fig |21 it is shown, how a network is build by ar- 
ranging the nodes. The minimum network forms a 
loop, which corresponds physically to a single magnetic 
bound state. A periodic continuation of this network 
leads to further adjacent loops, which get coupled by 
the nodes. In this way the network transmits poten- 
tials and the potential distribution is calculated by an 
appropriate iteration procedure. Currents are obtained 
only after obtaining the solution for the potential distri- 
bution within the whole network. Because of using just 
reflection and transmission coefficients so far, phase co- 
herence between adjacent loops is not explicitly included. 
Despite this fact, this network allows to simulate IQHE 
experiments with realistic sample geometry in almost all 
details [851 l9p| . The key-point of this approach is the 



representation of the back scattering function P, which 
will be considered in more detail in the main part of this 
paper. 

Considering the nodes of the network, they can be un- 
derstood as elementary QH - samples with a single EC 
pair experiencing back scattering, which is described by 
the function P. Originally this function had been used as 
representative for a QHE sample as a whole and it had to 
be defined for each LL separately as P — Pf Ai/)|68ll7l|. 
Av is the filling factor of the corresponding LL relative 
to half filling. The total sample behavior was obtained 
by summing up the components of the conductance ten- 
sor of all LLs. In 71] it had been shown, that in or- 
der to meet all symmetry relations found in experimental 
data [12, E3, P{Av) has to be an exponential function 
P(Av) — exp(—Av/k). On this background this func- 
tion has to be seen as a semi-empirical function so far. 
Nevertheless, just assuming the existence of this function 
allowed to explain already quite a number of universal 
features of the IQHE without referring to a special choice 
of the pre-factor k in the exponent. The function P(Av) 
also allows to obtain a scaling behavior by considering 
the factor k as a temperature dependent function k(T). 
In |7l| it has been shown, that any not necessarily known 
function k(T) is mapped out directly by the temperature 
scaling function for R xx and R xy - 

The basic idea, which subsequently led to the usage 
of P(Au) also within a network approach, is that cur- 
rent transport at high magnetic fields in the bulk re- 
gion may happen via tunneling between magnetic bound 
states. Such bound states are considered to be physically 
equivalent to ECs. Thus, these bound states are closed 
loops of directed channels and are created in real samples 
by smooth potential fluctuations at high magnetic fields. 
Tunneling between such loops happens preferably near 
the saddle points of the random potential. The most im- 
portant ingredient of our network is that the tunneling 
current can be handled as a back scattering process in the 
EC-picture like already outlined above. This makes our 
network approach esse ntial different from the Chalker- 
Coddington m odel I ll3l and a ll network approaches de- 
rived from this |115l Ill6t Ill7j . Just to clarify, the CC 
network maintains coherence on the entire network, while 
the largest coherent element in our network is a single 
loop, which at the same time defines the size of a net- 
work period L. If we consider the in-elastic scattering 
length Li n of the electron system and the size of the 
sample L p , then the CC network deals with the regime 
L < L p < L 

Lin l£ L < 

does not necessarily contradict the CC-network, but ad- 
dresses a different regime. As will be shown in this pa- 
per, the tunneling process is well re pres ented by the back 
scattering function P(Av) used in [83. It can be calcu- 
lated locally at the designated grid points from the mag- 
netic field and the local carrier density. On this basis 
it is possible to model the shape the sample by shap- 
ing the lateral carrier distribution. By doing so, it has 



L 



while our network considers the regime 
p . From this point of view our network 
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been shown for the first time, that the influence of sample 
geometry (such as contact arms) and carrier density in- 
homogeneities (like introduced by gate electrodes) can be 
successfully addressed on the basis of a network model in 
full agreement with experimental results |85l |90j . While 
the main task of |s3 was to present the technical aspects 
and examples for applications of the network model, the 
main purpose of this paper is to provide the theoretical 
and physical back ground. Since the network model is 
based on the idea of tunneling between magnetic bound 
states at saddles of the interior potential landscape, a sys- 
tematic study of this tunneling process will be presented. 
Also the interesting question, whether or not there is a 
correspondence between the edge channel picture and the 
bulk current picture in terms of mixed phases, will also 
be addressed. 



II. THE SADDLE POINT PROBLEM 

The problem of quantized transmission across a saddle- 
point constriction with magn etic field has been addressed 
in the past by Buttiker |114j . Biittikers paper was moti- 
vated by the theoretical discussion of quan tized conduc- 
tion steps discovered by van Wees et al[2(j and Wharam 
et al.[21| in split gate constrictions of a 2-dimensional 
electron gas. Since constrictions in these experiments 
are electrostatically induced with a pair of split gates, 
the potential is a smooth function (without hard walls) 
and the bottleneck of the constriction therefore forms 
a saddle. At this point our view of the role of poten- 
tial fluctuations meets the split gate situation referenced 
above. In this context we look for a link between the 
results of Biittiker for the bottleneck situation at high 
magnetic fields and our semi empirical back scattering 
function P(Az/). 

Following the paper of Biittiker for the high magnetic 
field limit, the trajectories of the states are given by 
equipotential lines defined by e = eV(x,y), where e de- 
notes the energy of the guiding center of the classical cy- 
clotron orbit, V the potential and e the electron charge. 
If Vq denotes the potential of the saddle, the case e < eVo 
describes a trajectory, which is repelled by the saddle, 
where the energy e is finally represented by the Fermi 
level. Using the considered saddle point as the center of 
our coordinate system, the minimum distance y n between 
trajectories of the same energy, which are classically re- 
pelled on opposite sides to the saddle point, is determined 
by e = eV(0,y n ) for e < eVo (see Fig. |3J). A similar sit- 
uation appears for e > eVb, but turned around by 90 
degrees and the minimum distance x n is determined by 
e = eV{x n ,Q). According to Biittiker, the transmission 
probabilities in terms of x n and y n in the high magnetic 
field limit for the case of e„ > is 




FIG. 3: Example for a two-dimensional potential modula- 
tion which forms saddle points. The minimum distances of 
the closed loops (magnetic bound states) are shown for two 
different cases: y n is achieved, if e is below the saddle energy 
eVo (bold line) and x n is achieved, if e is above the saddle 
energy. All possible loops represent equipotential lines and 
the selection of the active loops is made by the position of the 
Fermi level (e = Ef). 



and for the case e„ < is 



Tynn — 3jnn \ I T exp 



(3) 



In EqnEl and 03 the saddle energy has been defined 
as zero reference (Vo = 0). The underlying problem is 
basically a tunneling problem and S nm tells us that this 
tunneling process with the tunneling rate T mn happens 
mainly between the same states (states of same LL) on 
opposite sides of the saddle. The parameters u> y and uj x 
result from expanding the potential near the saddle point 
by 



V(x,y) « V + 0.5mw 2 x x 2 



0.5muj 2 y 2 



(4) 



Now we look for a relation between the transmission 
rates given above and the function P(Au), which repre- 
sents the ration between transmission and reflection R/T 
of edge states like used in the network model. If we look 
for x n while y n — (e„ = e(V — Vo) > 0), we can easily 
calculate x n from Eqn. 



^mn — fimn \ 1 ~f exp 



X t - 



(2) 



2e r , 



vriLOi 



(5) 
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FIG. 4: Contour plot of the artificial 2D Cosine potential. 
Dark indicates potential minima and the light colored areas 
indicates the potential maxima. L is the period of the poten- 
tial modulation and at the same time it is the grid period of 
the network. 



In the case e n < we get something analogous for y 2 : 



2e n 



(6) 



We introduce an artificial periodic potential modu- 
lation, which is a 2-dimensional Cosine - function and 
which has the same Taylor expansion like Eqn^l 



V(x, y) = V [cos{io y y) - cos(w x x)] 



(7) 



with V being the potential modulation. Fig. Q]shows a 
contour plot of the 2-dimensional Cosine potential modu- 
lation and we get a saddle point whenever the maximum 
of the Cosine in one direction meets a minimum of the 
Cosine in the other direction. In this way we get a peri- 
odic array of saddle-points. We consider one saddle to be 
located exactly at the origin of our coordinate system. If 
we think about a smoothly varying random fluctuation 
potential, the average curvature in both directions will be 
the same. In our representative potential in Eqn. we 
therefore simplify by setting 

(jj x — (jjy — (jj ^ 

which leads 

to uj x /ujy = 1. In order to re- write Eqn. [21 and Eqn. [3] we 
do some further substitutions. For the magnetic length 
Ib we use: 



I 2 - 



2neB 



(8) 



with h the Plank constant, e the electron charge and 
B the magnetic field. Using further the period length L 
of the periodic arrangement of saddle points 



we can rewrite Eqn. [5] and Eqn. |HJ 



L 2 e n 



!Jn 



2TT 2 eV 



T 2 F 

2n 2 eV 



(10) 



(11) 



Eqn. [S] - ^] now allow to re- write Eqn. [2] and Eqn. 
13 Since x n and y n are described exactly by the same 
function, but distinguished only by the case e n either 
above or below the saddle energy, we can rewrite Eqn. [21 
and Eqn. [3] into a single equation as follows: 



Pmn &mn \ 1 ~T" exp 



± 



L 2 e n eB 



eV h 



(12) 



The signes ± indicate the cases e n above or below the 
saddle energy. At this point it should be mentioned, 
that already Haug et al [25| have used a similar equation 
for addressing the problem of transport across magnetic 
bound states. They investigated experimentally the ef- 
fect of a gate electrode across the current path of a QHE 
sample at constant magnetic field. The major point of 
that experiment was to study the special situation of hav- 
ing integer filling outside the gate (pure EC transport) 
and sweeping the gate voltage through non-integer fill- 
ing below the gate. The authors describe the scattering 
through the bulk due to the presence of the gate by partly 
EC back scattering at the gate and partly back scattering 
via tunneling across magnetic bound states in the gate 
region. They achieved good agreement with the experi- 
ments already by a simplified network consisting of just 
2 or 3 magnetic bound states. 

In order to be able to compare with the formulation 
of edge channel back scattering in terms of P(Ais) , we 
have to realize, that the transmission process across the 
saddle T mn is formally a back scattering process in the 
EC picture. Therefore the transmission factor T mn in 
Eqn. EJcorresponds to a reflection R in the edge channel 
picture. We transform Eqn. ^Jinto P = R/T by using 
T mn as R and consequently using (1 — T mn ) as T. : 



P=* = 
T 



T 



1-T„ 



(13) 



Since the above tunneling process happens between 
edge states on opposite side of the saddle but the same 
Landau level (LL), d mn can be omitted for the further 
treatment. Further more we have to consider, that the 
energy e„ is determined by the Fermi level Ep. In this 
way we end up with P as an exponential function of the 
Fermi level Ep'- 



2tt/w 



(9) 



P = exp 



L 2 E F eB' 
eV h 



(14) 
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In this equation Ep is the Fermi level relative to the 
saddle energy. The ration L 2 /V can be understood as a 
measure of the " smoothness" of the potential modulation 
and eB/h is the well known number of states in a single 
LL per unit area. Since L is the distance between equiv- 
alent saddle points, L 2 is the area of the 2-dimensional 
basis cell of the network grid. It is interesting to note, 
that the pre- factor of Ep, which is (L 2 /eV)(eB /h), has 
the meaning of an average density of states (DOS) in the 
basis cell of the network, provided that we associate the 
LL broadening with eV . In order to consider this aspect 
in a little more detail, we ass ume a Gaussian shaped den- 
sity of states for a single LL |lOSj and calculate the DOS 
in the center of the LL as follows: 



DOS(e n = 0) 



eB 1 
h b\pH 



(15) 



where b is the energy broadening of the Gaussian peak. 
On this background we can calculate the change of the 
carrier density An in a single LL due to the change of 
Ep by AEp as follows: 



eB 1 
An = ——=AEp 

h by/ir 



(16) 



This is valid for a Fermi level close to the center of the 
LL, where the DOS is approximately constant. However, 
we can also represent An in terms of the filling factor by 
Av = An/riLL, where till = eB/h. This results in: 



Av = 



AE F 
~b~V^ 



(17) 



Using again the saddle energy as the zero reference, 
AEp becomes Ep and Av becomes the difference of the 
filling factor relative to half filling. Now we can substitute 
Ep by Av and rewrite Eqn. [21 



P = cxp 



-Av 



—^~-L — 
eV h 



(18) 



Since V is the amplitude of our artificial potential mod- 
ulation and b is the with of the representative DOS - peak, 
these two parameters correspond to each other. There- 
fore we can assume, that the pre-factors in the exponent 
will be a factor of the order of unity, just depending on 
details of the real potential fluctuations. On this basis 
we simplify 



eV 

and get the following exponential function: 



P = cxp 



-AvL 2 — 



(19) 



(20) 



It is easily seen, that the semi-empirical exponential 
function of Ref. [71| is completely recovered, except the 
fact, that we get a more complex pre-factor for the fill- 
ing factor Av. The explicit magnetic field dependence of 
the pre-factor is a new feature as compared to the earlier 
semi-empirical function. However, this has no impact on 
the major features of the function P(Av), just causing 
a magnetic field dependence of the width of the plateau 
transition regimes. If we take a closer look to the expo- 
nent, we find eB/h, which is the number of states per unit 
area in a LL and L 2 , which is the area of a grid period 
of the network. As a consequence we get a surprisingly 
universal result: The exponent of the function P(Av) 
represents the number of carriers being added (AN > 0) 
or removed (AiV < 0) from the half filled basis cell of 
the network grid. This happens if the Fermi level moves 
away from the center position of the LL. 



P = cxp [-AN] 



AN = AvL 2 — 



(21a) 



(21b) 



This is usually achieved by changing the magnetic field 
at constant carrier density. However, this change nor- 
mally happens continuously, meaning that the number of 
states in a single LL in one grid period is changed just in 
fractions of unity. This is possible, because in the regime 
close to a half filled LL there is a strong coupling between 
the individual grid periods so that the (single) electronic 
states are well de-localized and can be smeared over sev- 
eral grid periods. A change of the number of electronic 
states by AN = ±1 corresponds to a complete plateau to 
plateau transition. By considering the grid period length 
L as a measure of the smoothness of the potential fluctua- 
tions, we get a nice possibility for an interpretation of our 
result: A large L corresponds to a rather clean sample 
(high mobility) with very smooth potential fluctuations 
and a short L corresponds to a sample with strong po- 
tential fluctuations and hence, strong disorder (low mo- 
bility). As a consequence, a large L indicates a large 
area of the basis cell of the network and therefore only 
a small change of the filling factor Av and thus, only a 
small change of the magnetic field is needed in order to 
change the number of states by one. This means that 
we get sharp plateau to plateau transitions on the mag- 
netic field axis. In contrast, a small L, as expected from 
samples with strong disorder, leads to a small area of the 
basis cell, which then requires a large change of the filling 
factor Av for changing the number of states by the order 
of one. Therefore we get broad plateau to plateau tran- 
sitions in this case, exactly what we expect for samples 
with strong disorder. 

On the basis of geometric arguments and using a 
checkerboard as representative for the network (Fig. |5J) 
Eqn. 1211 allows also another simple interpretation: Us- 
ing directly the trajectories of the channels at different 
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FIG. 5: Spatial distribution of filled (dark colored) and empty 
(light colored) states for left: Ef > -Ell, middle: Ef = Ell 
and right: Ef < Ell- The dotted line indicates the basis cell 
of the network. 
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FIG. 6: Schematic representation of the effect of back scat- 
tering in the EC-picture. 



Fermi levels we get 3 cases as shown in Fig. (i) if 
the Fermi level is above the saddle energy (center of LL) , 
the filled states dominate the transport properties, be- 
cause they build an interconnected network, (ii) if the 
Fermi level coincides with the saddle energy, filled and 
empty states take up exactly half of the space, which 
also means exactly a half filled LL. (iii) if Ep is below the 
saddle energy, the empty states (insulating phase) build 
an interconnected network, leaving isolated droplets of 
filled states and leading to a high resistive regime. It 
should be mentioned, that already Mashida et al. used 
such a picture for the interpretation of their experimental 
resultsMIai- 



III. DISCUSSION 

By using a representative 2-dimensional Cosine poten- 
tial modulation, a result by Biittiker |114| for the saddle 
point problem at high magnetic fields was transformed 
into the back scattering function P(Av), which is used 
for the nodes of the network according to R ef. [8o| . One 
main extension as compared to Ref. |l!4j is the treat- 
ment of the tunneling current within the framework of 
the Landauer-Biittiker formalism. This means that a 
transverse current between a channel pair (tunneling cur- 
rent serves as a back scattering current between chan- 
nels), leads to a longitudinal potential drop in each of 
the channels. This can be considered as the essence of 
the EC-picture within the LB-formulation and there is 
no classical way to account for this behavior. In order to 
point out this fact somewhat more clearly, the situation 
is sketched in Fig. 




FIG. 7: Checkerboard model for bulk current transport in 
the QHE regime after Ruzin et al fl2| . The dark and light 
areas represent different quantum Hall liquids with different 
Hall conductance. Currents I and J are assumed to flow ex- 
clusively in the corresponding phase and the total behavior of 
the conductor is determined by the mixture of both phases. 



If there is a longitudinal current I L , it is carried 
dissipation-less by a channel pair and is represented on 
the right side by II = (pi — A 1 3)e 2 /h and on the left side 
by II = (P2 - ^i)e 2 /h, while (fii - /i 3 ) = (fi 2 - Hi) for 
current conservation reasons. A transverse current It, 
which may be enabled by some conduction or tunneling 
between opposite channels causes a longitudinal voltage 
drop. This is because It has to be represented also within 
the LB formalism according to It = (Ml — H2)s 2 /h on the 
upper edge and It = {^3 — Ha)g 2 /h at the lower edge, 
while (/ii — /12) = (/i3 — /U4) = A/j,l, again because of cur- 
rent conservation reasons. It is important to note, that 
any current (Il or It) within the EC-picture is defined 
only by a pair of directed channels, which must not be 
confused with a difference of opposite directed classical 
currents. As a consequence, it is not possible to assign 
a current to any single directed channel and the occur- 
rence of a perpendicular current It does not change the 
longitudinal current II- This is a fact, which must not 
be misinterpreted as a violation of Kirchhoff 's law of cur- 
rent conservation. Using finally the ration It/Il as an 
equivalent for R/T of the LB formalism, it is easily shown 
that the longitudinal voltage drop A/i^ is coupled to the 
Hall voltage like Afi L = (/zj — /i 2 ) = (fii — fi 3 )R/T, with 
R/T identified as P, which is the function derived in this 
paper. 

The above summarized view of the EC-picture has 
been used already in Ref. 60] in order to explain an 
anomalous behavior of the magneto transport of high 
mobility quasi 3-dimensional PbTe wide quantum wells. 
Also in Ref. [7l| this interpretation of the EC-picture 
has been successfully used which finally led to the ap- 
plication within the discussed network model. From this 
point of view the network model can be considered to be 
based purely on the EC-picture and the LB-formalism. It 
is therefore quite surprising, that also the ohmic regime 
is quite well represented. However, this puzzle can be 
resolved by reconsidering Fig. [S] which appears like a 
checker board representation of the network. From this 



point of view it meets the approach of Ruzin et al |l2j . 
Ruzin et al used such a checker board model for bulk 
transport in the QHE regime in order to obtain a univer- 
sal relation between longitudinal and transverse conduc- 
tances in the QHE and Fig, 5 summarizes this concept: 
Ruzin et al considered two competing quantum Hall liq- 
uids in the presence of a long-range random potential, 
with the magnetic field determining the area fractions 
close to 1/2. Fig. represents this situation for an arti- 
ficial periodic long-range potential. The "white" regions 
represent the phase with quantized Hall conductance ct\ , 
and the "dark" regions those of Hall conductance u-x- If 
the situation is such that the "white" phase percolates 
freely throughout the sample, the system will be on the 
plateau with a xy — <j\. Similarly if the "dark" regions 
percolate freely, the system will be on the plateau with 
&xy — Near the percolation threshold transport is 
controlled by quantum tunneling between clusters of the 
same phase at the saddle points. Ruzin et al related the 
net currents in the four quadrants to the four electric 
potentials at the edges V\ through V4 as given by: 

I = a 1 (V 2 -V 1 ) (22a) 



a = e /h a = 




a = 3e/h a = 2 e /h 




FIG. 8: Scheme for combining individual LLs, where each 
of them may have a mixture of quantum Hall liquid and in- 
sulating phase. In this case it is assumed that just the top 
LL is partly occupied and all lower ones are completely filled. 
The result is a combined system of mixed phases of different 
quantum Hall liquids, which means different Hall conductivi- 
ties. 



J = a 2 (V 3 -V 2 ) (22b) 



I' = <n(V 3 -V 4 ) (22c) 



J' = <r 2 {V i -V 1 ) (22d) 

Postulating current conservation within the individual 
phases (/ = /' and J = J'), Ruzin et al get automatically 
(y 2 — Vi) = (V3 — V4), which is identical to the results 
of our treatment as illustrated in Fig. [S] Further on, 
Ruzin et al obtained a universal semi-circle-relation be- 
tween longitudinal and transverse conductance. At this 
point it should be mentioned, that such a semicircle re- 
lation has also been obtained within our preceding work 
HH IZj| j which once more demonstrates the agreement of 
our results with the results of Ruzin et al. However, our 
model achieves additional results such as the filling fac- 
tor dependence and transport parameters in terms of R xx 
and R xy as well as the distribution of the applied voltage 
across the sample area. There is also a qualitative dif- 
ference between Ruzin's and our picture: In our picture 
the space is divided between occupied and non-occupied 
states, while Ruzin et al are talking about 2 competing 
quantum Hall liquids (QHL) of different Hall conductiv- 
ity. We can overcome this disagreement by recognizing, 
that in our model each involved LL is represented by 
a separate network with its own mixture of empty and 
filled states, which we associate with the insulator and a 
QHL phase of Hall conductance a — e 2 /h respectively. 



In the case of the IQHE there is always just one partly 
filled (top) LL and all lower ones are completely filled. 
For each LL the QHL phase exhibits a quantized Hall 
conductance of a = e 2 /h, while the Hall conductance 
for the insulator phase a = 0. For the complete sys- 
tem in the bulk current representation we have to sum 
up all contributions of all LLs locally, which leads to a 
situation like sketched schematically in Fig. [H] Putting 
forward a bulk current representation, one has to add 
all Hall conductivities of all involved LLs locally and the 
resulting picture becomes again similar to that one of 
Ruzin et al. At this point we want to make clear, that 
at this stage our network model does not include non- 
linear effects, which means that it considers the elec- 
tron system to be sufficiently close to thermal equilib- 
rium. However, for thermal equilibrium some similarity 
between bulk and edge current picture might not be that 
surprising. A recent theoretical treatment of non-linear 
tran sport through mesoscopic systems by Sanchez et al 
103J shows a possible magnetic field asymmetry, which is 
not within the scope of our network model in the present 
state of development. Nevertheless it is possible to de- 
scribe selective current injection into edge stripes also 
within our network model, like e.g. achieved by selective 
EC reflection due to gate electrodes [8^, while we cannot 
see how this could be achieved within the bulk current 
picture. Another important aspect to be mentioned at 
this point is that in our network the phase coherence is 
not maintained over more than one single grid period, 
but still exact quantization of the QHE is obtained. This 
is in agreement with theoretical results of F. Gagel et al 
[5o| who found a surprising stability of the QH plateaus 
against dissipation and phase destroying events. 
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FIG. 9: Sample layout as used for the numerical Simulations. 
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FIG. 10: Results for magneto transport data in the high 
field regime calculated with Eqn. 1201 as the coupling function 
for the nodes of the network. The calculation is done for 2 
different effective g-factors (g e ff = 0,5 and g e ff = 1.0) and 
two different values for L — 130nm and L = 185nm. 



IV. SIMULATION RESULTS 

If we use Eqn. | 20l as a replacement of the original 
function of Ref. [85|, we get results for the magneto 
transport simulation as shown in Fig. ^| and Fig. ^2 
The layout of the sample is shown in Fig. Consider- 
ing the individual plateau transitions at high magnetic 
fields, there is no qualitative difference to the results of 
the previous version. The major difference and improve- 
ment concerns the low magnetic field range. At very low 
magnetic fields the zeros in R xx disappear and the Hall 
effect starts as a straight line, and turns gradually over 
to a plateau behavior at higher fields. The reason for this 
is the explicit magnetic field dependence of the exponent 
in Eqn. I2UI which leads to wide plateau transitions at 
low magnetic fields and sharp plateau transitions at high 
magnetic fields. This results in an overlap of the ohmic 
(dissipative) regimes of several LLs at low fields. In order 
to get the individual filling factors in case of LL overlap, 
the individual spin split LLs have been calculated using 
the standard parameters for G aAs and a LL broadening 
according to T = To^/B |l08j | has been used by setting 
r = Q.bmeV. 

The actual calculation consists of two main parts: (I) 
For calculating the occupation numbers standard proce- 
dures are used and the Fermi level in the bulk region 
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FIG. 11: Results for magneto transport data in the low field 
regime calculated with Eqn. 1201 as the coupling function for 
the nodes of the network. The calculation is done for 2 dif- 
ferent effective g-factors (geff = 0, 5 and g e ff = 1.0) and two 
different values for L — 130nm and L — 185nm. 



is calculated by filling up the density of states (DOS) 
with the constant bulk carrier density. The DOS is com- 
posed by the superposition of the magnetic field depen- 
dent Gaussian shaped DOS of spin split LLs. The Fer- 
milevel for regions of non-zero electrostatic bare poten- 
tial (edges and gate regions) is forced to match the ob- 
tained (magnetic field dependent) Fermi level in the bulk. 
This allows to get a self-consistent electrostatic potential 
and a rearrangement of the carrier density at the edges, 
which results in an electrostatic edge potential accord- 
ing to Chlovskii et al [2|. (II)The such obtained lateral 
carrier density profile is transferred to the nodes of the 
network and in another self-consistent iteration proce- 
dure the lateral distribution of the bias voltage, which is 
introduced via the current contacts, is calculated. From 
this the potential difference for any designated pair of 
voltage probes is obtained as a function of the magnetic 
field. A metallic contact is realized by interconnecting all 
channels of all involved LLs at the position of the desig- 
nated contact point [g^- Extended metallic contacts are 
realized by creating arrays of such contact points. The 
current at the current contacts is calculated only after 
arriving at the self-consistent solution in the network, 
which allows finally to calculate the various resistances. 
This means, that in principle a constant supply voltage is 
used in the network model instead of a constant supply 
current like in most experiments. However, for calcu- 
lating the resistances for a standard QHE setup with a 
single current source this makes no difference. For achiev- 
ing a constant current mode, the potentials at the current 
contacts are additionally varied in a proper way during 
the iteration procedure in order to get the required pre- 
set current. In this way a one-to-one simulation of ex- 
periments using constant supply currents and measuring 
voltages at several potential probes is possible. 

Just for demonstration, the influence of spin splitting 
by using different effective g-factors such as g* = 0.5 and 
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g* = 1.0 and different values for L, such as L = 135nm 
and L = 185nm, is simulated. The results are shown in 
Fig. EHfor the high magnetic field range and in Fig. ITT1 
for low magnetic fields. For L — 135nm and g* = 0.5 spin 
splitting is hardly to observe even at high fields around 
5 Tesla, while for the case L = 185nm a well resolved 
spin splitting can be observed even for a small g-factor 
of g* = 0.5. At magnetic fields around 2.5 Tesla spin 
splitting is only weakly pronounced even for g* = 1 and 
L = 185mn. In summary one can say, that this over- 
all behavior looks quite realistic and a more systematic 
study on the possibility of a more accurate determination 
of an enhanced effective g-factor from transport data is 
planned separately. 



suits into a filling factor dependent back scattering func- 
tion P(Av). This function is used as the coupling func- 
tion for the nodes of a recently developed new network 
model for the QHE regime of 2D systems. As a result, 
the numerical simulation based on this network deliv- 
ers transport data, in which R xy starts with a classical 
straight line while R xx shows no zeros. With increas- 
ing magnetic field the IQHE sets-in with Shubnikov-de 
Haas oscillations in R xx , followed by zeros in R xx and 
plateaus in R xy , like seen in experimental data. For equi- 
librium conditions we have found an equivalence between 
the bulk current picture in terms of mixed phases like 
used by Ruzin et al and our network representation on 
the basis of the Landauer-Buttiker formalism. 



SUMMARY 



VI. ACKNOWLEDGEMENTS 



Starting with a treatment of the saddle point problem 
in the high magnetic field limit by Buttiker we have been 
able to show, that it is possible to transform Buttikers re- 



The author thanks M. Buttiker for stimulating discus- 
sions and comments. 



[23] 



M. Buttiker, Phys. Rev. B38, 9375 (1988) 

D.B. Chlovskii, B.I. Shklovskii, L.I. Glatzman, 

Phys.Rev. B46, 4026 (1992) [24; 

H. Hirai, S. Komiyama, Jpn. J. Appl. Phys. 34, 4321 

(1995) 

T. Christen, M. Buttiker, Phys. Rev. B53, 2064 (1996) [25 
H.P. Wei, D.C. Tsui, M.A. Paalanen, A.M.M. Pruisken, 
Phys.Rev.Lett 61, 1294 (1988) [26 
S. Koch, R.J. Haug, K.v. Klitzing, K. Ploog, Phys. Rev. 
Lett 67,883 (1991) [27 
S. Koch, R.J. Haug, K.v. Klitziung, K. Ploog, Phys. [28 
Rev. B 46, 1596 (1992) 

H. P. Wei, L.W. Engel, D.C. Tsui, Phys. Rev. B50, 14609 [29 
(1994) 

F.Hohls, U.Zeitler, R.J.Haug, R.Meisels, K.Dybko, 
F.Kuchar Phys.Rev.Lett. 89, 2768011-2768014 (2002) [30 
B. Huckestein, Rev. Mod. Phys. 67 (1995) 357 
S.L. Sondhy, S.M. Girvin, J. P. Carini, D. Shahar, Rev. [31 
Mod. Phys. 69, 315 (1997) 

I. Ruzin and S. Feng, Phys Rev, Lett 74, 154 (1995) [32; 
K.v. Klitzing, G. Dorda, M. Pepper, Phys. Rev. Lett. 
45, 494 (1980) [33 
F.F. Fang, P.J. Stiles, Phys. Rev. B27, 6497 (1983) 

F.F. Fang, P.J. Stiles, Phys. Rev. B 29, 3749 (1984) 

B.E. Kane, D.C. Tsui, G. Weimann, Phys. Rev. Lett. [34 

59, 1353 (1987) 

P. Streda, J. Kucera, A.H. MacDonald, Phys. Rev. Lett [35 
59, 1973 (1987) 

H.A. Fertig and B.I. Halperin, Phys. Rev. B 36, 7969 

(1987) [36 

J.K. Jain, S.A. Kivelson, Phys. rev. Lett 60, 1542 (1988) 

B.J. van Wees et al., Phys. Rev. Lett 60, 848 (1988) [37 

A. Wharam et al., J. Phys. C 21, L209 (1988) 
S. Komiyama, H. Hirai, S. Sasa, S. Hiyamizu, Phys. [38 
Rev. B 40, 12566 (1989) [39 ; 

B. J. van Wees, E.M.M. Willems, C.J.P.M. Harmans, 



C.W.J. Beenakker, H. van Houten, J.G. Williamson, 
C.T. Foxon, J.J. Harris, Phys. Rev. Lett 62, 1181 (1989) 

B. J. van Wees, L.P. Kouvenhoven, C.J.P.M. Harmans, 
J.G. Williamson, C.E. Timmering, M.E.I. Broekaart, 

C. T. Foxon, J.J. Harris, Phys. Rev. Lett 62, 2523 (1989) 
R.J. Haug, J. Kucera, P. Streda, K. von Klitzing, Phys. 
Rev. B 39, 10892 (1989) 

B. W. Alphenaar, PL. McEuen, R.G. Wheeler, R.N. 
Sacks, Phys. Rev. Lett. 64, 677 (1990) 

, T. Martin, S. Feng, Phys. rev. Lett 64, 1971 (1990) 

G. Muller, D. Weiss, S. Koch, K. v. Klitzing, H. Nickel, 
W. Schlapp, R. Losch, Phys. Rev. B 42, 7633 (1990) 

A. K. Geim, P.C. Main, P.H. Beton, P. Streda, L. Eaves, 

C. D.W. Wilkinson, S.P. Beaumont, Phys. Rev. Lett. 67, 
3014 (1991) 

J.A. Simmons, S.W. Hwang, D.C. Tsui, H.P. Wei, L.W. 
Engel, M. Shayegan, Phys. Rev. B 44, 12933 (1991) 

C. J.B. Ford, S. Washburn, R. Newbury, CM. Knoedler, 
J.M. Hong, Phys. Rev. B 43, 7339 (1991) 

P.C. van Son, F.W. de Vries, T.M. Klapwijk, Phys. Rev. 
B 43, 6764 (1991) 

H. Nii, M. Ohsawa, S. Komiyama, S. Fukatsu, Y. Shi- 
raki, R. Itoh, H. Toyoshima, Surface Science 263, 275 
(1992) 

S. Komiyama, H. Hirai, M. Ohsawa, Y. Matsuda, S. 
Sasa, T. Fujii, Phys. Rev. B 45, 11085 (1992) 

B. W. Alphenaar, A. A.M. Staring, H. van Houten, 
M.A. A. Mabesoone, C.T. Foxon, Phys. Rev. B 46, 7236 
(1992) 

S. Komiyama, H. Hirai, M. Ohsawa, Y. Matsuda, S. 
Sasa, T. Fujii, Phys. rev. B 45, 11085 (1992) 
R. Merz, F. Keilmann, R.J. Haug, K. Ploog, Phys. Rev. 
Lett. 70, 651 (1993) 

D. L. Maslov, D. Loss, Phys. Rev. Lett. 71, 4222 (1993) 
D.H. Lee, Z. Wang, S. Kivelson, Phys. Rev. Lett 26, 
4130 (1993) 



10 



[40] J.E. Muller, Phys. Rev. Lett. 72, 2616 (1994) 
[41] P.C. Main, A.K. Geim, H.A. Carmona, C.V. Brown, 
T.J. Foster, R. Taboryski, P.E. Lindelof, Phys. Rev. B 

50, 4450 (1994) 

[42] H. Hirai, S. Komiyama, Phys. rev. B 49, 14012 (1994) 
[43] A.J. Peck, S.J. Bending, J. Weis, R.J. Haug, K.v. Kl- 

itzing, K. Ploog, Phys. Rev. B 51, 4711 (1995) 
[44] K. Oto, S. Takaoka, K. Murase, Jpn. J. Appl. Phys. 34, 

4332 (1995) 

[45] R. Schuster, K.Ensslin, V. Dolgopolov, J. P. Kotthaus, 

Phys. Rev. B 52, 14699 (1995) 
[46] A. Kristensen, C.J. Kennedy, P.E. Lindelof, M. Persson, 

Semicond. Sci. Technol. 10, 1315 (1995) 
[47] B.J. van Wees, G.I. Meijer, J.J. Kuipers, T.M. Klapwijk, 

W. van de Graaf, G. Borghs, Phys. Rev. B 51, 7973 

(1995) 

[48] R.J.F. van Haren, W. de Lange, FA. P. Blom, J.H. 

Wolter, Phys. Rev. B 52, 5760 (1995) 
[49] I.V. Zozulenko, FA. Maao, E.H. Hauge, Phys. Rev. B 

51, 7058 (1995) 

[50] B.L. Johnson, A.S. Sachrajda, G. Kirzenow, Y. Feng, 

R.P. Taylor, L. Henning, J. Wang, P. Zawadzki, P.T. 

Colderidje, Phys. Rev. B 51, 7650 (1995) 
[51] T. Mashida, H. Hirai, S. Komiyama, Y. Shiraki, 

Phys.Rev. B54, 16860 (1996) 
[52] E. Yahel, D. Orgad, A. Palevski, H. Shtrikman, Phys. 

Rev. Lett 76, 2149 (1996) 
[53] M. Rahman, J.H. Davies, LA. Larkin, M.C. Holland, 

A.R. Long, J.G. Williamson, Phys. Rev. B 54, 16409 

(1996) 

[54] R.G. Mani, Journal Phys. Soc. Japan, 1751 (1996) 
[55] F. Gagel, K. Maschke, Phys. Rev. B 54, 13885 (1996) 
[56] J.Oswald, G. Heigl, M. Pippan, G. Span, T. Stellberger, 

D.K.Maude, J.C.Portal, Surface Science 361/362, 525 

(1996) 

[57] J.Oswald, G. Heigl, G.Span, A. Homer, P.Ganitzer, 
D.K.Maude, J.C.Portal, Physica B 227, 360 (1996) 

[58] Y.-C. Kao, D.-H. Lee, Phys. Rev. B 54, 16903 (1996) 

[59] T. Mashida, H. Hirai, S. Komiyama, Solid State Com- 
mun. 103, 441 (1997) 

[60] J.Oswald, G. Span, Semicond. Sci. Technol. 12, (1997) 
345 

[61] J. P. Bird, M. Stopa, K. Connolly, D.P Pivin, Jr., D.K. 
Ferry, Y. Aoyagi, T. Sugano, Phys. Rev. B 56, 7477 
(1997) 

[62] D. Shahar, D.C. Tsui, M. Shayegan, E. Shimshoni, S.L. 

Sondhi, Phys. Rev. Lett. 79, 479 (1997) 
[63] K. Tscmekhman, V. Tsemekhman, C. Wexler, D.J. 

Thouless, Solid State Commun. 101, 549 (1997) 
[64] R.G. Mani, Appl. Phys. Lett. 70, 2879 (1997) 
[65] J.Oswald, G. Heigl, G. Span, A. Homer, P. Ganitzer, 

D.K.Maude, J.C.Portal, Proc.of 12th Int. Conf. on the 

Appl. of High Magnetic Fields, Wrzburg 1996, Vol 1, 

227, World Scientific, Singapore (1997) 
[66] J.Oswald, G. Span, A. Homer, G. Heigl, P. Ganitzer, 

D.K. Maude, J.C. Portal, Solid State Commun. 102, 

391 (1997) 

[67] D. Shahar et al, Solid State Commun. 102 (12)(1997) 
817 

[68] J. Oswald, G. Span, F. Kuchar, Phys. Rev. B 58 (1998) 
15401 

[69] T. Mashida, H. Hirai, S. Komiyama, Y. Shiraki, Physica 

B 249-251, 128 (1998) 
[70] A. Furusaki, Physica B 249-251, 430 (1998) 



[71] J. Oswald, Physica E 3 (1998) 30 

[72] T. Mashida, H. Hirai, S. Komiyama, Y. Shiraki, Solid- 

State Electron. 42, 1155 (1998) 
[73] J.Oswald, G.Span, A. Homer, G. Heigl, P.Ganitzer, 

D.K.Maude, J.C.Portal, Proc. of the 8th Int. Conf. 

On Narrow Gap Semiconductors, Shanghai 1997, World 

Scientific, 304-307 (1998) 
[74] J.Oswald, G.Span, A. Homer, G. Heigl, P.Ganitzer, 

D.K.Maude, J.C.Portal, Proc. of the 8th Int. Conf. 

on Narrow Gap Semiconductors, Shanghai 1997, World 

Scientific, 300-303 (1998) 
[75] P.Ganitzer, A. Homer, M.Lucyshyn, B.Jamnig, 

J.Oswald, D.K.Maude J.C. Portal, Physica B 256- 

258, 600-603 (1998) 
[76] A. Homer, P.Ganitzer, G.Span, J.Oswald, Superlattices 

and Microstruct. 25, 191 (1999) 
[77] J. Oswald, A. Homer, P. Ganitzer, Microelectron. Eng. 

47 (1999) 31 

[78] T. Heinzel, Y. Acremann, K. Ensslin, E. Gini, H. Mel- 

chior, M. Holland, Physica E7(3-4), 804 (2000) 
[79] T. Mashida, S. Ishizuka, K. Muraki, Y.Hirayama, S. 

Komiyama, Physica E 6, 152 (2000) 
[80] M.T. Woodside, C. Vale, PL. McEuen, C. Kadow, 

K.D. Maranowski, A.C. Gossard, Phys.Rev. B 64 (4), 

041310/1-4 (2001) 
[81] H. Fujioka, S. Katsumoto, Y. Iye, Jap. J. Appl. Phys. 

40(3B), 2073 (2001) 
[82] K. Arai, K. Oto, S. Takaoka, K. Murase, Physica E. 

9(2). 243 (2001) 
[83] E. Ahlswede, P.Weitz, J. Weis, K. von-Klitzing, K. 

Eberl, Physica B 298(1-4), 562 (2001) 
[84] T. Mashida, S. Ishizuka, S. Komiyama, Koji Muraki, Y. 

Hirayama, Phys.Rev B63, 045318 (2001) 
[85] J. Oswald, A. Homer, Physica E 11, 310 (2001) 
[86] J. Oswald, Proc. 25th Int. Conf. on the Physics of 

Semiconductors in Osaka 2000, Springer Proceedings in 

Physics 87, 977 (2001). 
[87] T. Machida, S. Ishizuka, T. Yamazaki, S. Komiyama, 

K. Muraki, Y. Hirayama, Phys. Rev. B 65(23), 233304 

(2002) 

[88] E. Ahlswede, J. Weis, K. von-Klitzing, K. Eberl, Phys- 
ica E 12(1-4), 165 (2002) 

[89] A. Wurtz, R. Wildfeuer, A. Lorke, E.V. Deviatov, V.T. 
Dolgopolov, Phys.Rev.B 65, 075303(2002) 

[90] J.Oswald, 15th Int. Conf. on High Magnetic Fields in 
Semiconductor Physics, SEMIMAG, 5-9 Aug, 2002, Ox- 
ford, U.K., CD-ROM 

[91] J. Oswald ,Y. Ochiai, NAoki, L.-H. Lin, K. Ishibashi, 
Y. Aoyagi, Microelectronic Engineering 63, 91 (2002) 

[92] K. Ensslin, Superlattice and Microstruct 33, 425 (2003) 

[93] S. Komiyama O. Astafiev T. Machida, Physica E 20, 43 
(2003) 

[94] Y. Kawano, S. Komiyama, Phys. Rev. B 68, 085328 
(2003) 

[95] M. Blaauboer, Phys. Rev. B 68 (20), 205316 (2003) 
[96] K. Takashina, R.J. Nicholas, B. Kardynal, N.J. Mason, 

D.K. Maude, J.C. Portal, Phys Rev B 68 (23), 235303 

(2003) 

[97] M. Biittiker, P. Samuelsson, E.V. Sukhorukov, Physica 

E 20, 33 (2003) 
[98] E. Peled, D. Shahar, Y. Chen, E. Diez, D.L. Sivco, AY. 

Cho, Phys Rev Lett 91(23), 236802 (2003) 
[99] T. Machida, T. Yamazaki, K. Ikushima, S. Komiyama, 

Appl. Phys. Lett 82(3), 409 (2003) 



11 



[100] Y.Kawano, S. Komiyama, Phys. Rev. B 68(08), 085328 
(2003) 

[101] K. Arai, S. Hashimoto, K. Oto, K. Murase, Phys. Rev. 

B 68, 165347 (2003) 
[102] A. Cresti, G. Grosso, G.P. Parravicini, Phys. Rev. B 

69(23), 233313 (2004) 
[103] David Sanchez and Markus Biittiker, Phys. Rev. Lett. 

93, 106802 (2004) 
[104] P. Samuelsson, E.V. Sukhorukov, M. Biittiker, Phys. 

Rev. Lett. 92(2), 026805 (2004) 
[105] E. Peled, Y. Chen, E. Diez, D.C. Tsui, D. Shahar, D.L. 

Sivco, A.Y. Cho, Phys Rev B 69, 241305(R) (2004) 
[106] F.W. Hehl, Y.N. Obukhov, B. Rosenow, Phys Rev Lett 

93(9), 096804 (2004) 
[107] E.V. Deviatov, A. Wurtz, A. Lorke, M.Yu.Melnikov, 

V.T. Dolgopolov, D. Reuter, A.D. Wieck, Phys Rev B 

69(11), 115330 (2004) 
[108] M.Russ, A. Lorke, D. Reuter, P.Schafmeister, Physica 



E22, 506 (2004) 
[109] Y. Seo, B. Eom, I. Yu, K. Park, S. Lee, Solid State 

Commun 130(6), 391 (2004) 
[110] S. Kicin, T. Pioda, T. Ihn, K. Ensslin, D.D. Driscoll, 

A.C. Gossard, Physica E 21, 708, (2004) 
[111] R. Landauer, IBM J.Res.Dev. 1, 223 (1957) 
[112] M. Biittiker, Y. Imry, R. Landauer, S. Pinhas, Phys. 

Rev. B 31, 6207 (1985) 
[113] J.T.Chalker, P.D.Coddington, J. Phys. C 21 (1988) 

2665. 

[114] M. Biittiker, Phys. Rev. B 41, 7906 (1990) 
[115] S. Cho, M.P.A. Fisher, Phys. Rev. B 55, 1637 (1997) 
[116] D.P. Arovas, M. Janssen, B. Shapiro, Phys. Rev. B 56, 
4751 (1997) 

[117] P. Cain, R.A. R6mer, M. Schreiber, M.E. Raikh, Phys. 
Rev. B 64, 235326 (2001) 



